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Data from the Sudbury Neutrino Observatory have been used to constrain the lifetime for nucleon 
decay to "invisible" modes, such as n — > 3za The analysis was based on a search for 7-rays from the 
de-excitation of the residual nucleus that would result from the disappearance of either a proton or 
neutron from 16 O. A limit of tw > 2 x 10 29 years is obtained at 90% confidence for either neutron or 
proton decay modes. This is about an order of magnitude more stringent than previous constraints 
on invisible proton decay modes and 400 times more stringent than similar neutron modes. 

PACS numbers: 11.30.Fs, 11.30.Pb, 12.20.Fv, 13.30.Ce, 14.20.Dh, 29.40.Ka 
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Experimental signatures of the grand unification of 
the electroweak and strong interactions have been sought 
with increasing sensitivity for the past twenty-five years. 
Much effort has gone to identifying specific decay modes 
of free protons and bound nucleons as signatures of grand 
unification, but no signal has been observed to date. De- 
cay modes more unusual than those typically explored 



cannot, however, be ruled out (see, for example, [l|). A 
recent paper has even suggested a model in which n — > 3v 
becomes the dominant mode 0. Thus, the search for any 
mode which may have been missed by previous experi- 
ments is of fundamental interest. In this paper, the Sud- 
bury Neutrino Observatory (SNO) Q is used to search 
for what we will refer to as "invisible" decay channels, 
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i.e., those in which no visible energy is deposited in the 
detector via direct production of energetic, charged par- 
ticles. 

The search utilizes SNO's unique detection capabilities 
for low energy 7-rays, based on the Cherenkov light pro- 
duced by the resulting Compton-scattered electron. A 
generic signature for nucleon disappearance in 16 O arises 
from the subsequent de-excitation of the residual nucleus 
[4,5]. Approximately 45% of the time, the de-excitation 
of either 15 0* or 15 N* results in the production of a 7- 
ray of energy 6-7 MeV. SNO detects these 7-rays with 
good efficiency. In fact, the primary energy calibration 
source used by SNO is the 6.13-MeV 7-ray produced in 
the de-excitation of 16 N. 

A background to this nucleon-decay signal in SNO re- 
sults from neutral current (NC) interactions of solar neu- 
trinos. This is due to the 7-rays of similar energies that 
are produced as a result of neutron captures on nuclei in 
the detector. In Phase I (D2O), neutrons were detected 
through observation of the single 6.25-MeV 7-ray result- 
ing from neutron capture on deuterium. This capture 
efficiency is 0.29, but the threshold applied to the anal- 
ysis to limit low energy backgrounds reduces the over- 
all neutron detection efficiency to e„ = 0.144 ± 0.005 
6]. In Phase II (D 2 0+NaCl), two tonnes of NaCl were 
added to the one kilotonne of D2O. Neutron captures on 
35 C1 release 8.6 MeV of energy in 7-rays, with most cap- 
ture events producing multiple 7-rays. The correspond- 
ing capture efficiency in Phase II is 0.90 and, due to a 
relatively high analysis threshold, the overall neutron de- 
tection efficiency is e' n = 0.399 ± 0.010 7]. The multiple 
7-rays from neutron captures in Phase II result in a more 
isotropic distribution of Cherenkov light, which can be 
used as a further discriminant for identifying the neutron- 
induced component. However, the principal advantage in 
comparing Phase I and Phase II data lies in the fact that 
7-rays from the nucleon-decay signal are detected with 
similar efficiencies in SNO, while neutrons produced by 
8 B solar neutrinos are detected with very different effi- 
ciencies. These characteristics are used in what follows 
to measure an upper limit for nucleon disappearance. 

In terms of the data from Phase I of SNO, the rate 
of nuclear 7-ray production can be related to the ap- 
parent production rate of neutrons by taking account of 
the detection efficiencies for neutrons, 7-rays and particle 
misidentification as follows: 

i? 7 e 7 / 7 „ = R n — £ n 7Z-NC 

where i? 7 is the rate of nuclear 7-ray production due 
to nucleon decay; e 7 is the efficiency for detecting the 
nuclear 7-rays above the analysis energy threshold; f jn 
is the fraction of the detected nuclear 7-rays which are 
mistaken for neutrons; R n is the extracted neutron de- 
tection rate nominally attributed to NC interactions; e n 
is the neutron detection efficiency for the fiducial volume 
and analysis energy threshold; and TZnc is the actual 



production rate of neutrons due to solar neutrino NC 
interactions. Similarly, for Phase II data, 
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where AR n is the difference between the extracted neu- 
tron detection rate attributed to NC interactions in 
Phase I and that implied by data from Phase II. 

In order to compare Phase I and Phase II rates under 
the same assumption for the underlying CC spectrum, 
results from SNO data were used in which the CC com- 
ponent was constrained to follow the shape of a standard 
B energy spectrum Table I summarizes the rele- 
vant results from these two phases. The extracted num- 
bers of CC, NC and ES (elastic scattering) events include 
the subtraction of all known backgrounds (as detailed in 
and 0), including atmospheric neutrino interactions 
which might identically mimic the nucleon-decay signal 
via the removal of a proton or neutron from 16 O. 

The number of NC events extracted in Phase I [(| 
was 576.5 for a livetime of 306.4 days, yielding a rate 
of R n = 686.8 ± 83.9 per year, with statistical and sys- 
tematic uncertainties added in quadrature. Similarly, for 
Phase II, 1265.8 NC events were implied based on a live- 
time of 254.2 days, yielding a rate of R' n = 1817.6 ±136.6 
per year. Thus, accounting for the relative neutron de- 
tection efficiencies, 



ARn 



rt n 1 ^n 

(686.8 ±83.9) 
30.8 ±97.3 



(656.0 ±49.3) 



Thus, an upper limit of AR n < 180.6 per year at 90% 
confidence limit is obtained using a standard, Bayesian 
prescription (which is also in good agreement with fre- 
quentist prescriptions) [fj. 

From 4], a vanishing neutron from the 16 O nucleus 
results in an excited state which has a branching ratio 
of 44% for producing a 6.18-MeV 7 and 2% for a 7.03- 
MeV 7. For a vanishing proton, the distribution is nearly 
the same, with a branching ratio of 41% for a 6.32-MeV 
7 and 4% for a 7.0-MeV 7. The signal extraction pro- 
cedures previously used for solar neutrino analyses were 
applied to simulated nuclear 7-ray lines of these ener- 
gies, combined with a simulated solar neutrino signal. 
The numbers of additional NC events extracted relative 
to the actual NC signals generated were then expressed 
as fractions of the generated nuclear 7-ray signals. The 
values of / 7 „ and /' were then determined by the ap- 
propriate weighting of these fractions in accordance with 
the relative branching ratios given above. For Phase I 
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TABLE I: Signal extraction results for CC constrained to B shape. Error bars are the quadrature sum of statistical and 
systematic uncertainties. 



Analysis Parameter 


Phase 1 (pure D2O) 


Til TT* / T\ /~\ 1 TVT /"II \ 

Phase 11 (D2O + NaCl) 


Fiducial Volume 


6.97 x 10 8 cm 3 


6.97 x 10 8 cm 3 


Energy Threshold 


T eff > 5 MeV 


T cfi > 5.5 MeV 


Livetime 


306.4 days 


254.2 days 


CC Events 


f967.7±117.9 


1430.3±97.1 


ES Events 


263.6±29.2 


163.7±23.8 


NC Events 


576.5±70.4 


1265.8±95.2 


Neutron Detection Efficiency 


0.144±0.005 


0.399±0.010 


NC Event Rate (R n & R' n ) 


686.8±83.9 yr" 1 


1817.6±136.6 yr' 1 


Equivalent Phase I NC Rate (Rn & ^R' n ) 


686.8±83.9 yr" 1 


656.0±49.3 yr" 1 



* The Phase II data set used for this analysis is identical to that presented in 7] 



data, it was found that f yn = 0.99+Q02 for both neutron 
and proton decay modes. This is as expected since the 
neutron signal in pure D2O results from a 6.25-MeV 7- 
ray, which is virtually indistinguishable from cither 6.18 
MeV or 6.32 MeV within the energy resolution of the 
detector. The distributions are, therefore, nearly 100% 
covariant. For Phase II data, /' = 0.75^0;^ (again, 
nearly identical for either decay mode). Once more, this 
is roughly what is expected given the additional isotropy 
information. The lower value of /' derived reflects a 
compromise within the fitting procedure between provid- 
ing a good description of the isotropy distribution and 
the energy spectrum expected for neutrons. These same 
simulated nuclear 7-ray lines were also used to determine 
e 7 and e 7 . For neutron (proton) decay modes, these were 
found to be 0.51 ± 0.01 (0.59 ± 0.01) and 0.361 ± 0.005 
(0.425 ± 0.006), respectively. 

Thus, an upper limit can be deduced for the number 
of decay 7-rays at greater than 90% confidence level of 
Rh m < 443 per year for neutron decay and i?Jj m < 385 
per year for proton decay. An upper bound to invisible 
modes of nucleon decay can now be established as follows: 



where N np is the number of neutrons or protons (depend- 
ing on decay mode) within the D2O fiducial volume which 
are bound in 16 (1.85 x 10 32 ), and e 7 is the efficiency 
for the decay to result in the release of a 6 or 7-MeV 7- 
ray (0.46 for neutron modes and 0.45 for proton modes). 
Therefore, the comparison of Phase I and Phase II data 
from SNO implies that, at greater than 90% confidence 
level, 

for neutron modes: Tj nv > 1.9 x 10 29 years. 



Prior to this paper, the best constraint on n —> 3v used 
by the Particle Data Group was based on Kamiokande 
data in which higher energy, but much weaker, branches 
of the de-excitation of the oxygen nucleus were consid- 
ered and yielded a limit of r > 5 x 10 26 years |]|. It has 
been proposed that a similar analysis could be carried out 
with data from Super- Kamiokande and, by making use 
of the carbon nucleus, possibly even in the KamLAND 
detector 0. It has also been noted that the disappear- 
ance of a proton from the deuteron in heavy water de- 
tectors would result in a free neutron, which could then 
be captured to yield a detectable signal for invisible pro- 
ton decay (see, for example 0]). This has already been 
used to yield a lower bound to the proton lifetime in 
excess of 10 28 years for such modes [6,12]. Lead perchlo- 
rate has also been suggested as a possible future detec- 
tor medium to search for invisible nucleon decay, making 
use of de-excitation of the nuclear hole that would be 
left in 35 C1, with an estimated sensitivity on the order of 
10 30 years for a one kilotonne detector [13,14]. Owing to 
the extremely low background levels in SNO, the prin- 
cipal branches of the de-excitations for 16 O have been 
probed here and have yielded limits which are within a 
factor of 5 of this level. Thus, the constraint presented 
here is about an order of magnitude more stringent than 
the recently published limits on invisible proton-decay 
and 400 times more stringent than previous limits on 
neutron modes, such asm "iv. 
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for proton modes: r- mv > 2.1 x 10 years. 
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